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Superconductivity of quasi-one-dimensional organic conductors with a quarter-filled band is investi- 
gated using the two-loop renormalization group approach to the extended Hubbard model for which both 
the single electron hopping tj_ and the repulsive interaction Vj_ perpendicular to the chains are included. 
For a four-patches Fermi surface with deviations to perfect nesting, we calculate the response functions for 
the dominant fluctuations and possible superconducting states. By increasing V_\_, it is shown that a d-wave 
(singlet) to /-wave (triplet) superconducting state crossover occurs, and is followed by a vanishing spin 
gap. Furthermore, we study the influence of a magnetic field through the Zeeman coupling, from which a 
triplet superconducting state is found to emerge. 
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1. Introduction 

Superconductivity in quasi-one-dimensional (quasi- ID) or- 
ganic conductor, (TMTSF^X, has been studied extensively 
in the conditions where charge and spin fluctuations play an 
important role due to the low dimensionality of the Fermi 
surface. 1 ' 2 ' The possibility for triplet state superconductiv- 
ity, besides the singlet one, is an issue of current interest 
in these materials. 3 " 6 ' From recent NMR measurements on 
(TMTSF)2C104, it has been suggested that spin-triplet super- 
conductivity may emerge out of a singlet state under mag- 
netic field. This can occur as a field-induced phase transition 
that can compete with a FFLO state under strong magnetic 
field. 5 ' It is important for such a study to take into account the 
influence of low-dimensional fluctuations due to the strongly 
anisotropic band structure of (TMTSF^X. The existence of 
spin fluctuations is supported by the NMR experiments in the 
normal phase of these materials. 5,7 ' 

Several theoretical works have been devoted to the field- 
induced phase transition to the triplet superconducting (SC) 
state. Shimahara 8 ' 9 ' pointed out such a transition using the 
pairing interactions mediated by spin fluctuations, which 
have both components for the spin-singlet and triplet pair- 
ings. Combining mean-field and RPA methods, Belmechri et 
a l 10, li) nave snown (h a t un( j er fjeid a singlet-triplet super- 
conducting transition is possible, along with the occurrence 
of a FFLO state at intermediate strength of the magnetic field. 
A similar transition has been shown to occur by Aizawa et 
a/., 12 ' 13 ' using the RPA method for the extended Hubbard 
model that includes intersite repulsive interaction for longitu- 
dinal and transverse directions along the chains, and the Zee- 
man coupling of spins to a finite magnetic field. As shown for 
the Hubbard model with on-site repulsive interaction, the RPA 
approach, 14 " 16 ' which sums up a higher order of perturbation 
for electron interactions, suggests the importance of the pair- 
ing interactions mediated by spin fluctuations. Thus it is of 
interest to further examine fluctuations of both density waves 
and superconducting pairings, where the spin gap is essential 
to the existence of a singlet superconducting state. 

These features can be properly taken into account by the 



renormalization group (RG) method. 1718 ' The effect of a 
magnetic field on low-dimensional systems has been studied 
by the RG method mainly for the one-dimensional cases. 19 ' 20 ' 
Noticeable progress has been achieved in studying supercon- 
ductivity in zero field for the case of quasi- ID systems with 
many chains. 21 " 23 ' However, the different mechanisms by 
which triplet superconductivity can be stabilized in such sys- 
tems, especially in finite magnetic field, have not been fully 
investigated within the RG scheme. 

In the present work, we use the RG method up to the two- 
loop level 24 ' to study the competition between the li-wave sin- 
glet SC state (SCd) and the /-wave triplet SC state (SC/) in 
quasi- ID systems with interchain electron hopping and repul- 
sive interactions. It is demonstrated that superconductivity is 
driven by the interplay of interchain interaction and the nest- 
ing deviations. The crossover from the SCd state to the SC/ 
state occurs with increasing the interchain repulsive interac- 
tion in the presence of nesting deviations. The effect that a 
Zeeman coupling to the magnetic field can have on the sta- 
bility of the SGi state and the emergence of a triplet SC/ 
state is also studied in details. In §2, we give the formula- 
tion of the RG technique for a many-chains quasi- ID sys- 
tem at quarter-filling, in the presence of magnetic field and 
nesting deviations. Using a four-patches decomposition of the 
Fermi surface, we derive the flow equations for the SCd and 
SC/ response functions in the superconductivity channel, and 
for the spin-density-wave (SDW) and charge-density-wave 
(CDW) responses in the staggered density-wave channel. In 
§3, the results for the possible states as a function of inter- 
chain Coulomb interaction, nesting deviations and magnetic 
field are presented. The conditions for the stability of the SC/ 
triplet state in the calculated phase diagrams are given. Sum- 
mary and discussion are presented in §4. 

2. Formulation 

2.1 Model 

In order to study the superconductivity for the (TMTSF^X 
salt, we consider the quasi- ID extended Hubbard model, 



2 J. Phys. Soc. Jpn. 



Full Paper 



K. Kajiwara et al. 



(a) 



fb) 



Fig. I. (Color online) (a) Extended Hubbard model with intra- and in- 
terchain repulsive interactions. Interchain hoppings are shown by dashed 
lines, (b) The quasi-lD Fermi surface. The locations of the patches, 4 in 
number, are shown by the closed circles. 



given by 



where 
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Here Cj^ a [c: l a ) is the annihilation (creation) operator of an 
electron at the site j and chain /, with spin (7 ((7 =t,4X and 
n j.i.o = c ) i c iJ,o- As illustrated in Fig. 1 (a), the quantities 
t\\,t±\,t±2,h,ll,U and V± are in order the nearest-neighbor 
intrachain and interchain hoppings, next-to-nearest-neighbor 
interchain hopping, Zeeman field = 1), chemical poten- 
tial, and finally the on-site and interchain repulsive interac- 
tions. 

Equation (1) is rewritten by making use of the Fourier 
transform, c jjM = {IJ*x)~ ll2 'L h e ik \\ j+m c a (k), where the 
lattice constant is taken as unity. The kinetic energy, eq. (2) 
can then be written as Ho = Lfc.o- £(&)co-(fc)c a (fc). Here N± is 
the number of chains, L is their length, k — (k\\ , k) and e(k) = 
— 2f||Cosfc|| — 2t±i cosk — 2fx2Cos(£|| + k). In the quasi-lD 
case, we have the open Fermi surface since | , |fj_2 1 <C fii. 
The Fermi surface, which is a function of k, is divided into two 
parts for right-going (p = +) and left-going electrons (p = —) 
[see Fig. 1 (b)]. Further by adopting the linear-dispersion re- 
lation, eq. (2) is rewritten as 

H = £ e Pi a(k)cl a (k)c Pta (k), 

k,a 



(4) 



e,,a(fc) = v[p*|i -*£„(*)], 



u h 2t\ 2?2 

F n (k) = k p + a— H cosk — p — sinA:, 

2v v v 



(5) 
(6) 



where p = +/- (or R/L), k° = n/A, v = V2f|| , a = +/- (=t 

/ 4), fi = fxi +t±.z/ V2, and ?2 = txz/vZ- Terms of <9(? 2 ) and 
k dependence of the Fermi velocity are discarded. 



Expressing the interaction in terms of forward and back- 
ward scattering, 17 ' eq. (3) is rewritten as 

2kv ^ 



Hi = 



'{a}(q,k u k 2 ) 



- 1 - ki,k2,q,{o} 

x 4,(7! (fel)cL,cr 2 (fel -g)4,a 3 ( fc 2 -q)cR,a 4 (k 2 ), (7) 

where 
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The couplings G 2 ± and G|| stand for the forward scattering 
with spin being anti-parallel and parallel, respectively. The 
amplitude Gi± denotes that of the backward scattering with 
anti-parallel spins. The Sfcp-Umklapp scattering due to the 
quarter-filling is discarded. The coupling constants depend on 
the wave vector perpendicular to the chains, 21 ' and the defini- 
tion is the same as in ref. 24. The bare scattering amplitudes, 
which correspond to the initial values for the RG equations, 
are given by 



1 



{U + 2Vlcosq), 



r 21(qM,h) 



2nv 

— [U + 2V[cos(k l -k 2 )]., 



2nv 



= — [2Vlcos(Jki-Jk 2 )-2V^cos9]. 



(9a) 
(9b) 
(9c) 



In the following, we shall only retain the backscattering part 
(VY) and neglect the forward scattering contribution Vj_ of 
eqs. (9a)-(9c); VT is known to be involved in the enhancement 
of 2k° CDW fluctuations. 23 - 25 ' 

2.2 RG equations for the vertex couplings 

We consider the partition function, which is represented in 
the path integral form, 



S[¥*M 



(10) 



where S[\f/*,\j/] — Sq[w*,W] +Si[W*,Y] is the action corre- 
sponding to the Hamiltonian (1). The fields y/-*' are the 
Grassmann variables for the electron degrees of freedom 
and Dy/*^\j/ corresponds to the integration measure for the 
Grassmann variables. In the Fourier-Matsubara space, the free 
and interacting parts of the action So[\j/*, y/] and Si[y*,y] are 
respectively given by 18 ' 

s [y*,y]= £ [4, a (MM~V;U(fc)vWfc)> C 11 ) 

fc, iO) n , (7, p 



Si[Y*,V] 



H?[Y*,Y], 



(12) 



where T is the temperature (kg = 1 throughout this work) and 
//j G [y/*, ip] is obtained by substituting the fermion operators 
for the Grassmann variable e.g., c p a (k — q)—> V/ p ,o-(fe — q), 
etc., where k = (ku,k,CH,) and q = (q^q,(O m ), with (0 n and 
C0 m being the Matsubara frequencies for fermions and bosons, 
respectively. The Green's function for the free fermions is 
given by g° p CT (fe, ico„) = [i(0„ — e p , a (k)] . In order to exam- 
ine the behavior at low temperature, we proceed to implement 
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Fig. 2. (Color online) Peierls (a) and Cooper (b) bubbles. Solid (broken) 
line represent right (left) going electron, p = +(— ); the square is the in- 
teraction vertex i.e., G[±, G 2 ±, and Gii. Slashed lines represent electrons 
on high-energy shell Ei+m < \e„ a (k)\ < £";, while crossed lines represent 
electrons at higher energy or lower / values. 



the successive partial integrations of eq. (10) on high-energy 
shells. This leads to the renormalization of the inner shell ac- 
tion (5<) at step I of the procedure, where / > is a RG param- 
eter for the scaled energy E\ = Ee (2E[ is the scaled band- 
width at Z). 21 ) The partial integration of high-energy outer 
shell degrees of freedom is performed perturbatively with re- 
spect to So. The use of the linked cluster expansion allows one 
to write the result in the form 



£^y/*^y/exp 



(13) 



where (• • • )o. c denotes the contribution of connected diagrams 
to the outer energy shell of width E/dl, where dl <C 1. The 
action S[ + di, at the step l + dl, thus contains additional renor- 
malization terms, which for the coupling constants read 



& G v{q,k x ,k 2 ) (0 — G v(q,ki ,k 2 ) (l + dl)- G V (q,k u k2) W 



0Lr v(qM,k2) 



5G "v(qM,k 2 ) 



(14) 



where v — 1 _L,2 _L, and ||. The logarithmic contributions at 
the one-loop level (n = 2) for SG^? 2 ^ ^ are given by the 
diagrams shown in Fig. 2. 

This renormalization consists of two parts, which come 
from the Peierls channel and Cooper channels, namely 



0U v(q,k h k 2 ) 



HqM,k 2 ) 

The Peierls bubble Fp (q, ki,ki) is given by 
2nvT 



t C(qM,k 2 )- 



(15) 



P(q,ki,k 2 



LNi 



LLL 



-L Vi,V 2 k i(0„ 



where 
•A 



V,Vi,V2 
P(q,k,ki,k 2 ) 



V,Vi ,Vo 

(—l) m p G 



Vi(q,ki,k)( J v 2 {q,k,k 1 )i 



(16) 



(17) 



and the summation of Vi, V 2 (— 1 -1,2 _L, ||) is taken for the 
fixed v, as shown explicitly later. The quantity m p ' l ' 2 de- 
notes the number of the permutation of Grassmann variable. 
££ s ' is performed in the outer shell region. The summation 
is written as 

£ =2£©(|e L)CT ,(fc-< z )|-^), (18) 

k k 

where YSt' denotes the summation in the region of Ei + di < 
\£R,a(k)\ < E/ and 0(x) is the Heaviside step func- 
tion with the definition ©(0) = 1/2. In the RG proce- 
dure, the external momentum is fixed at the Fermi sur- 



face, i.e., (Ap CT] (/q),/q), (-k^ a ,(ki - q)M - q) for the in- 
coming states and (k§ a2 (k 2 ),k 2 ), (-lc^^(ki-q),ki-q) 
for the outgoing states. Thus for the momentum summa- 
tion of the bubble with q = (q^,q), the longitudinal mo- 
mentum qn is determined by the momentum conservation 
for the respective vertex (Vi,V2). 24) Using (LN^y 1 ^™' = 
(2nvN±)~ l Y.k So s dER^ a (k), and performing the Matsubara- 
frequency summation, eq. (16) is rewritten as 



P Hq,hk) 



— Y Y r v ' v i'' 



-L Vi ,v 2 k 



v 2 

kuh) 



de R}0 (k) 



f(£R,a (k)) - f(e L , a > (k-q)) 



£ Rt a(k)-£ L , a i(k-q) 
x&(\e L>a ,(k-q)\-Ei), (19) 

-El+dl 



where J os de = de + JI^" de. The Fermi 



distribution 

function/^) in eq. (19) is given by f(x) = l/[exp(x/T) + 1], 
which will be treated in the low-temperature limit. Following 
the RG procedure one can rescale the energy 2E/ up to the 
original band width 2E, so that all the energies A(= t\,t 2 ,h) 
are rescaled according to 



A(l) = Ae' 



(20) 



where the anomalous corrections to A due to the two-loop 
corrections are small and have been neglected for the present 
choice of parameters (see Appendix). At zero temperature, eq. 
(19) is calculated as 

1 



r P(qMM) 



LV J v <Vl,V 2 
L J P{q,k,ki,k 2 ) 



-L Vi,v 2 k 

x 2 L J P(q,k,ki, 



\dl, 



where 



and 



2E 



iPiq.k.kj.o-Oj) 



^ + \Y^ ua _ ai (i)y 



(21) 



(22) 



q,k,kj,0~Oi 



(i) 



■e Rt a(k)-e Li(J ,(k-q) 
= + 2t\ (I) [cos k + cos(k — q)] 

— 2t\ (/)[cos^,- + cos (A;; — q)\ 

— 2?2(/)[sinA; — sin (k — q)] 

+ 2t 2 (l) [sin A:; + sin (kj — q)] 

+h(l)(a-Gi). (23) 

In the above equations, a =t and I are replaced by o = + 1 
and —1, respectively. In a similar way, the Cooper bubble is 
calculated as follows: 



r C(q,k u k 2 ) 



2%vT 



LN 



LLL 

-L V],V2 k i(On 



X J, 



V,V[,V2 



C(q,k,ki,k 2 ) 

where q' = q — k± — k? and 

V,V] ,v 2 



g% a (k,ico n )g La ,(-k - q' ',-i(On), (24) 



V,Vl,V2 

C(q,k,k u k 2 ) 



(-l) m P 



G V[ (q-k 2 +k,ki ,k) G v 2 (q-k\ +k,k,k 2 ) ' 



(25) 



4 J. Phys. Soc. Jpn. 



Full Paper 



K. KAJIWARA et al. 




" I I " ~ "4 i" 

Fig. 3. (Color online) The RG equations for the coupling constants at the 
one-loop level. The slashed and crossed internal lines represent the elec- 
trons on the high-energy shell and within the shell, respectively. The di- 
agrams with the slash and cross lines interchanged (not shown) are also 
taken into account. 



and the summation L£/ s '" is given by substituting £l,c'{~^ ~ 
q') for e i CT /(fc — q) in eq. (18). The Fermi surface of the 
external momentum is given by (k^ ffi (k\ ) , k\ ), (—kp a , (k.2 — 

q),k2 — q) for the incoming states, and (^5 (^2)^2)1 
(— kp 1 (k\ —q),k\ — q) for the outgoing states. The longitudi- 
nal momentum vector q 1 ^ is also determined by the momentum 
conservation for respective vertex. Performing the Matsubara- 
frequency summation, eq. (24) is rewritten as 

2 



FX 



* Vi,V 2 



JV -L k Vi " 

x / de R , a (k) — — 

Jos. eR,a{k) + e Lia i{-k-q ') 

x&(\e L ^(-k-q')\-E,). (26) 

At zero temperature, eq. (26) is calculated as 

1 



r C{q,h,k 2 ) ~ N U L, J C(q,k,h,k 2 ) 
X Vl,V 2 k 

X 9 H Ic(q',k.kj,o-Oj}> 
L i=l,2 



where 



and 



Ic(g , ,k,k i ,a-a i ) 



IE 



(27) 



(28) 



Y ^k,k,.a-a,( l ) = e«,a(fe)+eL,a'(-fe-90 

= +2fi(/)[cosfc — cos(k + q')] 

—2t\ (I) [coski — cos(fc; + q')] 

— 2?2(Z)[sinA: — sin (k + q')} 

+2t 2 (l) [sin A:; - sin (Jfcj + 5')] 

+/j(/)(a-a ( -). (29) 

The functions h( g ,k,ki,o-Oi) and k{ q ,k,ki,o-Oi)> whi ch are equal 
to unity in the one-dimensional case and in the absence of the 
magnetic field, are reduced by t\ and ?2 in the quasi-lD case. 

The RG flow equations at the one-loop level are shown di- 
agrammatically in Fig. 3. Including the two-loop corrections 



derived in the Appendix, the RG equations take the form 



dl ^v{q,k u k 2 ) ~ 2fJ, L*~v(q,k,klte) 



1 _ , . . , »- _( 2 ) 



-^rG v {qMM)Y. z . 



J- q't' 

where S^ fcjti ^ is the contribution from the one-loop RG, 

(2) (2) 
Whlle 4(U,W,*0 ^ Z v(q,k U k 2 ,q>,k>) ^ tW °- l0 °P C011 - 

tributions coming from the self-energy and the vertex part, re- 



spectively (see Appendix). 24 ^ The quantity E 
ten as 

w ll(«,«li 2 ) — Lri - L (9^1>*:) Lr IIW.* : >*2) J P(^,* 1 ,jfc2) 

1) 



V{q,k,k^ki) 



is wnt- 



+ G ||( 9 ,A: 1 ,A:)Glx(< /A i :2 )/p (9ii . jl( . 2ii . i) 



l-L(?-*2+^l,S;) G 2_L( ? -fei+*:Afe)''c( ? -*:i-*2A*i,*2) 
2_L {q-k 2 +k,k i ,k) lJ l±{q-k i +k,k,k 2 ) 1 C(q-k i -k 2 ,k,k 2 ,k\ ) ' 

(31a) 



— G 



C '2i.($,*,jfcl,fe) ( -'2±{q,k l ,k)G2±(q,k,k2) 1 p( (l ^k l ,k 2 ) 

r r(2) 
J 1_L {q-k 2 +kM ,k)^l±(q-ki +k,k,k 2 ) I C(q~k l ~k 2 ,k,h ,k 2 ) 



^2Uq-k 2 +kM .k)G2±(q-k l +k,k.k 2 Yc(q-k i -k 2 .k,ki,k 2 )' 

(31b) 

r r r(2) 

~||(3,(fc,*l,*2) _ U ' 1 - L (9 , ^l>*) U ' 1 - L (9A*2)''P(g^,* 1 ,jfc2) 
(0) 



r(0) 



+ G ||(?,ii,i) G ||(?M2) / P( ? Mt,fe) 



G || (q-h+k.ki ,k) G \\ {q-h +k,k,k 2 ) I C(q-k l -k 2 ,k,h MY (3 



where for A = P, C, 



/ (0) _ V r 

1 l(q,k,k 1 ,k 2 ) ~ 2 



(32a) 



('=1,2 



r (l) 



X(q,k,k u k 2 ) 4 



7 H [h(q,k,k 1 ,2r)+h(q,k,k 2 fl)]> ( 32b ) 



=±1 



L h(q,k,ki.2r) 



X(q,k,k u k 2 ) 4 



(32c) 



r=±li=l,2 



- 1 



Note that in the absence of interchain couplings, E^,/ ^ t , 
coincides with the expressions already obtained by Montam- 
baux ef a/., 20 - 1 at /; 7^ in the ID case. 

2.3 RG equations for response functions 

Now we calculate the response functions for CDW, SDW, 
SCd, and SC /. The composite fields of corresponding order 
parameters are defined as 



CDW \ 



(33a) 



fe.cr 
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(a) 




Fig. 4. (Color online) Renormalization group equations for the response 
functions (a), and order parameter vertex part for the source field (b). 



y SDWl 



(<?p) = \IYn7 L c R,o( k ) c L--o{k - q P ), 



(33b) 



-L k.a 



^scw(gc) 



LN 



^(<Tcos^)c ffi<7 (fc)c L ,_ (J (-fc + qc), 



-L k.a 



(33c) 



^SC/(<Zc) = \ljrr ^(sinA:||Cos^)c«, a (fc)c L . a (-fc + qt c ), 

" V -L k.a 

(33d) 

where the density-wave modulation qp = (2kp, it) is the nest- 
ing vector and qc = (0,0). 

For simplicity, we use the four patches model for the Fermi 
surface, for which k = 0, it, as shown in Fig. 1 (b). We thus 
end up with 12 independent coupling constants. Such a choice 
allows a qualitative description of the SC states, i.e., eqs. (33c) 
and (33d) having respectively the gap functions cosfc and 
sin A: | cos£. This corresponds to nodes located at k = ±it/2 for 
SCd, and ku = and k = ±it/2 for SC/. Hereafter, we use 
or it for the momentum perpendicular to the chain. While the 
?2-terms in eqs. (23) and (29) vanish for q,k,kj = Q,it, nest- 
ing deviations can be incorporated in the Peierls channel by 
introducing the following factor for the Peierls channel in eq. 
(30), 

E 

(34) 



L* = — 



-4(1) 



where t-£) is a cut-off energy. It follows that I t * ~ 1 for 
E > ?2 (/) and I,* ~ for E -C f| (/)■ A cutoff procedure sim- 
ilar to eq. (34) has been used in studying the density-wave 
problem in many-coupled-chains case. 28) 

We calculate the response functions by adding a linear cou- 
pling of order parameters to source fields in the action, 18 ' that 
is S/,[v, yf*,h,h*] — X^fyx^u +c.c), where hp is a source 
field in the /l=CDW, SDW, SCd, and SC/ channels. The total 
action at step I is given by 

S[w,W*Ah*] l =S[w*Mi 

where is a renormalization factor for the order parameter 
vertex (z^(0) = 1), and %u(l) is the response function with 
Xn(0) = 0- 

Up to the one-loop level (see Fig. 4), the RG equations for 
these quantities are 

^Zcdw (<?p) = ^cdw (<?p)4/2% , (36a) 



^ZSDw(gp) = — ZSDW(9P)%) 



dl 
d 
dl 



Xscd{qc) = —4cd(lc)h/2, 



^xscfiqc) = i:,4cfiic), 



TZV 
1 

1Z\ 



(36b) 
(36c) 
(36d) 



where the respective three-point vertices z^ obey the RG 
equations: 

d 



dl 



lnzcDw(gp) - [—Gix(k,o,o) - g i±(ti fi.ji) 



d 

dl 
d 
dl 



+ G\\ (jcfl,7c) + G\\ (7c,o,o)] 4/2% > (37a) 

lnzsDw(qp) =[G 2 x( 7r ,o, 7 r) + G2±(7[fi.o)]lq, (37b) 

m 2sc<i(<?c) =[-Gix(o,o,o) + G : ix(^,o,^) 

+ G 2±(7r,0,7r) _ G f 2X(0,0,0)]V 2 ' ( - 37c - ) 



^■lnzsc/(9c) =G||( W)0)W ) ~ G||( 0) o,i 



())• 



(37d) 



The initial values are Zfi |/=o = 1. In the above equations, If, is 
defined by 

E 



h = 



(38) 



E + hQ)' 

Note that eq. (37b) represents the flow equation for the 
transverse SDW. The longitudinal SDW response is ob- 
tained by substituting [Gj^^o) + Gi±( w ,o jW ) + G ||(w,0,w) + 
g ||(ji,o.o)]4/24] f° r me r.h.s. of eq. (37b). These two flow 
equations become equivalent in the absence of h. 

3. Singlet versus triplet superconductivity 

For the numerical calculations that follow, parameters of 
the model are fixed at U = 4t\\ and t\ = 0.2t\\ unless stated 
explicitly. The unit of the energy is fii, which is set to unity 
f|| =1, and we take E= 2t^ . The scaling parameter I is replaced 
by / = ]n(E/T), where T can be squared with the actual tem- 
perature introduced in §2. 

In the following study of the four-patches Fermi surface, 
we will also examine the spin gap A a , which is governed by 
the combination of coupling constants 29 ' 



G(j+ — 



1 



T 2± (0,0,0) 



— G 



(0,0,0) ■ 



|(jr,0,0). 



(39) 



For a non zero A<j, the coupling G a + takes a positive value 
at / = 0, but moves to a fixed point with a large negative 
value. At a qualitative level for the spin gap, we shall use 
A CT = E[ a (= Ee~ la ) where l a is determined by the condition 
G a+ (l a ) = -0.7. 24) Thus a non zero spin gap is obtained for 



the couplings, G 



2_L(0,0,0) 



<0,G 



(0,0,0) 



>0, G 



2±(ji,0.0) 



< 0, and 



G||(^ o o) > 0. Here we note that in the absence of magnetic 
field, the combination of couplings for the spin gap, eq. (39), 
can be rewritten as 

1 



' [g 



1±(0,0,0) • 



r l±(jr,0,0)J 



(40) 



From G a +, the uniform susceptibility within RPA takes the 
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Fig. 5. (Color online) Temperature dependence of response function with 
f| = 0.0, h = 0.0 for V\ = 0.0 (a), and 1.0 (b). Response functions X\i are 
normalized by noninteracting response function Xo(= l/icv). 




Fig. 6. (Color online) Possible states in the V^-T plane with f| = 0.0 and 
h = 0.0. The solid line represents the boundary of the respective phase, and 
the dotted line is the energy, t\ , corresponding to interchain hopping. The 
dashed line denotes the spin gap A obtained from eq. (39). 



form 30 - 31 ' 



Xs = 



XsO 



1 - 2jivG a+ Xso 



(41) 



where XsO denotes the susceptibility for free electrons given 
by27rvXo-tanh(£/2r). 

3.1 Case of perfect nesting (t% = 0) 

Since the RG equations in the case of perfect nesting re- 
duce to those obtained for two-coupled chains, 29 ^ the domi- 
nant state in the absence of magnetic field is then either SGi 
or CDW. For V] 3 = 0, the system is characterized by a singlet 
SCd state as shown in Fig. 5 (a). Such a result is well known 
for the Hubbard ladder model with repulsive interactions. 32 ' 
The SDW correlations, which are the most dominant fluctua- 




Fig. 7. (Color online) Phase diagram in the V\-h plane with h = 0.0002. 
The solid lines designate the boundary for respective phases. The line be- 
tween SDW and SCf refers to V!>„. 



tions at high temperature, becomes less dominant at low tem- 
perature due to the formation of a spin gap. For larger 
than a critical value, V"j\., the CDW becomes in turn the dom- 
inant state as illustrated in Fig. 5 (b). The SCf correlations 
are enhanced, but show a weaker increase than CDW at low 
temperature due to the spin gap. 

In Fig. 6, the most dominant states are shown in the V±-T 
plane. At high temperature (T > t\), corresponding to one- 
dimensional regime, SDW is dominant. With decreasing tem- 
perature, the effect of interchain hoppings grows and the dom- 
inant state is in turn either SCd or CDW. On a tempera- 
ture scale, this occurs at a characteristic temperature given 
by the solid line. This scale decreases and becomes zero at 
v l = v ±c(— 0-83); it increases monotonically for V] 3 > V± c . 
This behavior resembles to that of the spin gap, A CT , as shown 
by the dashed line in Fig. 6. The existence of such a critical 
point of Vj" c has been also shown for two-coupled chains. 29) 
A similar behavior shown by A a indicates that spin degrees 
of freedom are also critical at , The behavior around the 
quantum critical point, V£ [ s ascribed to the competition be- 
tween SCd and CDW, which are associated to different spin 
gaps. The spin gap for SCd is formed by interchain pairing, 
whereas that of CDW results from intrachain interactions. 
Therefore A CT vanishes at the critical point where the sym- 
metry of the gap changes. 

We now turn to the effect of magnetic field, h, which brings 
new states due to its influence on the spin gap. When A<j is de- 
stroyed by the magnetic field, the magnetic state is expected to 
be either the transverse SDW or SCf. Such a region actually 



exists when Vj 3 



The phase diagram in the V] 3 - 1\ plane 



is shown in Fig. 7. The boundary between SCd and SDW (and 
also between SCf and CDW) is estimated from the condition, 
A(j — > 0. It is found that the effect of h on SCd is larger than 
it is for CDW The SDW state thus takes place for V\ < Vh, 
whereas the SCf state appears for V] 3 > v£ . Moreover, v£ 
is a critical value for the crossover between SDW and SCf as 
the dominant fluctuations. In our two-loop approach, the fact 
that both SDW and SCf of Fig. 7 vanish for h = is at vari- 
ance with the emergence of SDW state found in ref. 25 at low 
temperature using one-loop RG. 

3.2 Nesting deviations and superconductivity 

We now examine the SC state in the presence of nesting 
deviations, which is the main subject of the present paper. 
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We look at the possible states at tX 7^ and low temperature 
by choosing V\ = 0.7 (a), 0.88 (b), and 1.0 (c), namely for 
small, intermediate, and large interchain couplings. Figures 8 
(a), (b), and (c) show the temperature dependence of coupling 
constants which give rise to the response functions for SCd, 
SCf, and CDW respectively. The other couplings (not shown 
in the figures) only contribute a lesser degree to the response 
functions. These functions are traced in Figs. 9 (a)-(c). 

For V\_ = 0.7 [Fig. 8 (a)], the values of G 
C?2_L(ff,0,ff) at the fixed points are positive, while those of 
C?1J_(0,0,0) and G2x(o,o,o) are negative. The SCd response func- 
tion is then the strongest and becomes the dominant state. 
Since the coupling constant G2X (0,0,0) changes its sign and 
becomes relevant, a spin gap appears [see eq. (39)], which is 
crucial to SCd. As for the relevant coupling G^x^o^) (> 0), 
it enhances both SDW and SCd, as seen from eqs. (37b) 
and (37c). Figure 9 (a) shows the temperature dependence 
of the response functions, where the SDW state is found to 
be the dominant state at high temperature, but becomes sub- 
dominant at low temperature. The amplitude of SDW corre- 
lations is reduced by nesting deviations. Thus it is found that 
SCd pairing is induced by spin fluctuations with the coupling 

G2±(k,0,k)(> 0). 



For Vj 3 = 0.88 [Fig. 8 (b)], the relevant couplings are given 
by G||(jj ,0 .ji)(> 0) and G||(o,o,o)(< ^) which are quite different 
from those obtained in (a) for weaker v£; SCf fluctuations are 
then dominant, as seen from eq. (37d). The spin gap vanishes 
since the combination of couplings given by eq. (39) remains 
positive. Thus if A a = 0, the dominant contribution to the SC 
coupling comes from density fluctuations for which the rel- 
evant coupling with parallel spins, namely Gii( W) o )W )(> 0), is 
connected to charge fluctuations. It is worth noting that long 
wave length spin fluctuations also promote the SCf state since 
that G11 (000) < strengthens SCf correlations and also the 
uniform spin susceptibility [eq. (41)], as it will be discussed 
later. The CDW fluctuations are well developed at tempera- 
tures just above the region where SCf is dominant [Fig. 9 (b)]. 
At these temperatures, the CDW response function is nearly 
the same as that of SDW, implying a coexistence of spin and 
charge fluctuations for intermediate Vj 3 . We finally note that 
the dominance of SCf state is obtained in the presence of 
nesting deviations (tZ 7^ 0), which suppresses the divergence 
of CDW in the low temperature limit. 

For yj 3 = 1.0 [Fig. 8 (c)], the relevant couplings 
are G UjiM)) (> 0), G| K ^ ^)(> 0), G l±{j[fifi) (< 0), and 
Gi±(jcfi,n)(< 0)> which lead to a CDW state as seen from eq. 
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4x10" 



Fig. 10. (Color online) The T - V* "phase diagram" with if/f] = 0.001 
and h = 0.0. 



(37a). The behavior of Gm( Wj o jW )(> 0) is similar to that of (b), 
and then a similarity of CDW and SC/ response functions is 
expected. However, the nature of the spin gap as induced by 
Gijfjf o,0) (< 0) differs. For large Vj\ it is expected that the 
dominant state is CDW, whereas SC/ is sub-dominant. Here 
the difference in sign for Gi jf^Q^) in comparison to the case 
(a) is essential to the predominance of either SGi or CDW. 

The dominant states are summarized in the Vj_-T phase di- 
agram of Fig. 10. This phase diagram shows some similar- 
ity with the one of Fig. 6 in the sense that the SDW state is 
found at high temperature, namely for T > f| . However, for 
T < f|> the SC/ state comes in between the regions of the SGi 
and CDW states. In the zone where T < 1 1 and > V± c , the 
CDW state is dominant and the SC/ state is sub-dominant, 
but the CDW state is suppressed and the SC/ state becomes 
dominant for f| ^ — a result well explained by the fact that 
nesting deviations have a detrimental influence primarily on 
density-wave correlations. The effect of t%, when < , is 
small, because the SCd state is less affected by nesting devia- 
tions in that region. 

At this point we would like to comment on the behavior of 
the uniform magnetic susceptibility, Xs> given by eq. (41). In 
Fig. 11, the temperature dependence of this quantity is shown 
for the Vl = 0.7, 0.88, and 1 .0 cases considered above, which 
correspond to the SGi, SC/, and CDW states, respectively. It 
should be noticed that % s increases when the system is enter- 
ing in the SC/ state, indicating an enhancement of long wave 
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Fig. 1 1 . (Color online) Temperature dependence of the uniform suscepti- 
bility for V\ = 0.7, 0.88, and 1.0, where the parameters are the same as 
Figs. 9 and 8. The arrow denotes the temperature at which the SC state 
becomes dominant with decreasing temperature. 
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Fig. 12. (Color online) Phase diagram of the most dominant states in the 
T-h plane with = 0.78 and fj/f] = 0.005. In the inset, the temperature 
dependence of the SCd (solid line) and SC/ (dotted line) response func- 
tions is shown where the thin (bold) solid line corresponds to the case for 
h = (0.0002). 



length spin correlations. This contrasts with the cases where 
SCd and CDW states prevail and % s decreases rapidly due to 
the formation of a spin gap. 

Now we study the effect of a magnetic field on the SGi 
state for the intermediate interchain interaction, namely in the 
region where the sub-dominant SC/ state is close to SGi. Fig- 
ure 12 shows the T vs h phase diagram for the most dominant 
states at yj 3 = 0.78. At low temperature, there is a crossover 
from SGi to the SC/ state due to the suppression of the spin 
gap with h. In the inset, the temperature dependence of the re- 
sponse functions for SCd and SC/ are shown at finite (h= 
0.0002) and zero magnetic field. The SCd correlations are 
strongly suppressed by h, while the magnetic field has essen- 
tially no influence on SC/ correlations. 

4. Summary and Discussion 

We have examined by the two-loop RG method the singlet 
SGi and the triplet SC/ superconducting states in the frame- 
work of a 4-patch model with nesting deviations. The triplet 
SC/ state, which is absent in the model for intrachain inter- 
actions only, is found to develop from the combined effect of 
interchain repulsive interactions and nesting deviations. 

It is of importance to return to the mechanism of formation 
of the superconducting state and examine the scattering pro- 
cesses that are pertinent to this state. This is shown in Fig. 13, 
where the symbol ± in the kn-k plane denotes the sign of the 
superconducting gap function for the SGi [i.e., cos A: in eq. 
(33c) (a)] and SC/ [i.e., sinfc|| cosfe in eq. (33c) (b)] phases. 
The scattering processes for parallel and anti-parallel spins 
are described by a dashed arrow, while the long continuous ar- 
rows stand for the nesting vector. The sign of the renormalized 
coupling constants are also stated in each case considered. A 
sign change of the gap function following the scattering oc- 
curs for positive renormalized coupling constants, while the 
sign remains the same for attractive ones. Thus all the scatter- 
ing processes shown in Fig. 13 gain the energy. 

Figure 13 (a) depicts scattering processes with anti -parallel 
spins leading to the SGi state. The interaction with interchain 
momentum transfer n, i.e., the interband scatterings, give rise 
to the repulsive couplings among which G2x(7r,o,7r) lS the dom- 
inant contribution and G 1± ^ %) a l so becomes relevant, lead- 
ing to the growth of spin fluctuations as explained after eq. 
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(38). The sign of the gap function changes in these scattering 
process. For the couplings G2x(o.o.o) ar, d ^ij_(o,0,0)> which be- 
come attractive through the renormalization, the sign of the 
gap function remains the same in the scattering process. 

We now turn to the SC/ state, where the relation between 
the sign change of the gap function and the forward scatter- 
ing process with parallel spins is shown in Fig. 13 (b) (upper 
panel). For the repulsive G^ n ^ n y the gap function is opposite 
in sign, whereas the sign is the same for the attractive G||(o,o,0)- 
The backward scattering shown in the lower (b) panel also fa- 
vors the SC/ state, since the signs of the respective coupling 
constants, —G^ k q^ < and — G|| (o.o.o) > 0, are consistent 
with the change of the sign in the gap function. 

We now comment on the crossover from the SCd to SC/ 
states as a function of V]\ The dependence of relevant cou- 
pling constants at low temperature (T = 10~ 3 ) is shown in 
Figs. 14 (a) and (b), which correspond to the SCd and SC/ 
cases, respectively. For the SCii state, the coupling Gu_(o,o,o) 
evolves from negative to positive values, suggesting a van- 
ishing spin gap. The resultant reduction of the SCd state by 
v£ is reasonable because the interchain spin singlet state is 
destroyed by the formation of CDW state for moderate Vj 3 . 
With increasing Vj 3 , the repulsive interactions G2±( K .o tJ[ ) and 
Gix(jr,o.7r) decrease and the amplitude of the SCd correlation 
is reduced. 

In Fig. 14 (b), the coupling Gii^o^) increases from zero, 
while the coupling G\\ (o,o,o) decreases to negative value. Nega- 
tive G|| (o,o,o) indicates the absence of spin gap. These features 
yield in turn the development of the SC / state, as a conse- 
quence of charge fluctuations and nesting deviations. From 
Fig. 14 (b), the Vj 3 interval, where G a + > , suggests that the 
spin gap vanishes for V] 3 ^ yt> However, it is found that the 
SC/ state moves to the CDW state by noting that G a+ > for 
V\ = 0.88, and G a+ < for Vj 3 = 0.92. 

Finally, we comment on the effect of Vj 3 on the SDW state. 
Within the conventional treatment of RPA, 14_16) the interfer- 
ence effect between the scattering at different transverse mo- 
menta is neglected. For example, in eq. (9b), which is rel- 
evant to the SDW state, the onsite-repulsion U is retained 
but the effect of V] 3 is strongly reduced due to the summa- 
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Fig. 14. (Color online) The dependence of the renormalized coupling 
constants at T = 10~ 3 for the SCd (a) and SC/ (b) states, and for G a +- 
Here f| = O.OOlfi and h = 0.0. The initial values for Gijmo o)> ''UfjiO,*)' 
G 2 ±(o,o,0). G 2± (s,0,s). G ||(0.o,0) and %jt,0,s)' me respectively given by 
0.45, 0.45, 0.45, 0.45, 0.00 and 0.00 for V\ = and 0.72, 0.18, 0.45, 0.45, 
-0.27 and 0.27 for = 1.2. The arrow represents V*. 



tion of &i = ±7t and k2 = ±7T. However, the present RG 
method shows a clear effect of Vj 3 on SDW as illustrated by 
the strength of density-wave correlations in the Vj 3 — U plane 
at fixed T (Fig. 15). The SDW and CDW regions are defined 
by the corresponding response functions that become larger 
than the bare value Xo at I = by a factor greater than 10 6 . The 
domain that separates the two regions shows a much reduced 
amplitude of the density-wave response functions, while its 
area reduces by the decrease of temperature. The region for 
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Fig. 15. (Color online) The regions of SDW and CDW dominant correla- 
tions in the Vj-U plane defined by the condition ^sdw(Xcdw) > 10 6 ^o at 
the fixed temperatures T = 2 x 10 (dashed line), 4 x 10~ 4 (dash-dotted) 
and 8 x 1(T 4 (solid line). Here *f = 2 X 1(T 4 . In the remaining in between 
region ^sdw(Zcdw) < 10 6 ^o- 



SDW is suppressed with increasing V± due to the increase of 
the charge fluctuations. 
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Appendix: Two-loop level renormalization group 

The renormalization of the interchain hopping t\ is given 

by 24 > 



i 



(A-l) 



-L q,k,K 

where the second term of r.h.s. comes from interactions. Such 
an effect is negligibly small for the quarter-filled band al- 
though the strong effect is expected due to the Umklapp scat- 
tering for the half-filled band, e.g., the confinement of the in- 
terchain hopping. 26 ' 27 ^ 

We show the results at the two-loop level in eq. (30). 
Applying the calculations already obtained at perfect nest- 
ing and half-filling 24 ^ to the case of the nesting deviation 
at quarter-filling, the self-energy and vertex corrections are 

4 2 (Uw.*'r s ?(L,w,*r whichread 



,(2) 



+ G T. n (c,',k 2 ,k') J Uq'M,k') + G l n (q'-k l +q,k>) J Hq'-ki+q,k') 
+ G Z n (q',-k 1 +q,k') J nq l -k 2 +q,k l )i (A-2a) 
^( 2 ) - J 

^l±(q,k u k 2 ,q',k') ~ J 2(q+k';k l ,k 2 ;k',k'-q') 

X [ G lMq+q',k l) k 2 )G2_L(q-k 2 +k',k'-q',k')G\\( q -k l +k' ,k',k'-q') 
+ GlX(q+tf M,k 2 ) G \\(q-k 2 +k' ,k'-q' ,k')G2_L(q-k l +k' ,k' ,k' -q')] 
+ J2(-k';-k l ,-k 2 ;n-k' ,Ji-k'+q') 
x [GzLfa-k! ,k l ,k l -q , )G\\{k 2 -k> \k 2 -q\k 2 )G\±( q - q i ,k l -q l ,k 2 -q') 



+ G 



\{k x -k\k x .k x -q')G2\_(k 2 -k' .k 2 -q' .k 2 )G\L{q-q' .k x -q> .k 2 -q'^ 



r(2) 



(A-2b) 

a 2i(q,ki,k2,q',k') = J 2(q+k l ;k u k 2 ;k',k'-q') 

x [G2±(q+q',k l ,k 2 )G2_L(q-k 2 +k',k'-q',k')G2±{q-k [ +k',k',k'-q') 
+ G2i_(i 1 +q' : k x ,k 2 )G\\(q-k 2 +k' ,k' -q' ,k')G\\(q-ki+k' ,k' ,k' -q') 
+ G\\(q + q> ) k l) k 2 )Gi_L(q-k 2 +k',k'-q',k')Gi_L(q-k l +k',k',k'-q')] 
+ J2(-k';-k l -k 2 ;it-k',Ti-k'+q') 

x [G2±(k l -k l ,ki,k i -q') G 2±(k 2 -k' .k 2 -q' ,k 2 ) G 2X{q-q' ,k\-q' .k 2 -q 
+ G\\(k l -k l \k x .k x -q')G\\(k 2 -k! .k 2 -q! ' .k 2 )G2L(q-q< ' .k x -q> '.k 2 -q') 
+ G\±(k l -k' .k l :k\-q')G\X{k 2 -k' .k 2 -q> ,kj)G\\( q -q' .k,-q' .k 2 -q')] 

(A-2c) 



r(2) 



& \\(q,h,h,q'^) = J n<l+k h M,k 2 -k>,k'-q>) 

x [G\\(q+qi,k l ,k 2 )G2±(q-k 2 +k',k'-q',k')G2±(c,-k l +k',k',k'-q') 
+ G \\{q+q',k i ,k 2 ) G \\(q-k 2 +k',k'-q',k')G\\(q-k- l +k',k',k'-q') 
+ G2_L(q+q>,k l ,k 2 )Gu_(q-k 2 +k',k'-q'.k')Gl_L(q-k l +k',k',k'-q')] 
+ J2(-k';-k l .-k 2 -n-k' ,n-k' +q') 

x [G2±(ki-ki ,k l ,k l -q')G2_L(k 2 -k' ,k 2 -q' ,k 2 )G\\(q-q' ti-q' .k 2 -q<) 

+ G\\ fa -k',ki ,k { -q')G\\ [k 2 -k! .k 2 -q' .k 2 )G\\ (q-q' ,k x -q'.k 2 -q') 

+ G l ^ kl _ k ik i k i -q')G\]_(k 2 -k',k 2 -q l ,k 2 )G2i_(ci-q'.k i -q l ,k 2 -q l 

(A-2d) 



where G l n (q, k ,k') is s ivenb y 



G ln(qAk>) = E G l(9,k^) 



(A-3) 



v=l±,2±, 



™&h{q,k,vy J \{qfrk>) ^ g iven as follows. For \Y* kk , (l)\ < 
E, 



J, 



0(q,k,k') 



'l(q,k,k') 



2E\n[ 



q,k,k'S 

4£ -^,oO- 



16E 2 



(A-4a) 
(A4b) 



For|yP (Z)|>£, 



T 1rl r 4E + \ Y q,k,k'A l )\l , v p / 1\\ / a c \ 

Jo(qM') = 2Eln I TB i ivP 7Kl\ s % n ( Y qAk>fi( 1 )), ( A ' 5a ) 



./ 



2E 



IE 



\{q,k,k>) 



(A-5b) 



The quantity J2{ q +k"-M,k 2 \k' ,k") is also S iven b Y 

J2{q+k";k l .k 2 ;k',k") = j [j \(q+k" -k u k' -k") + J\{q+k"-k 2 ,k' ,k")] ■ 

(A-6) 

In the present calculations carried out at the two-loop level, 
we treated nesting deviations as the dominant effect. The in- 
fluence of magnetic field at that level remains to be examined. 
This will be the subject of a separate publication. 
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